The Morphology of the Galactic Dark Matter Synchrotron 
Emission with Self-Consistent Cosmic Ray Diffusion Models 

Tim Linden/'Q Stefano Profumo/'^'Q and Brandon Anderson^'^ 

^Department of Physics, University of California, 
1156 High St., Santa Cruz, CA 95064 
Santa Cruz Institute for Particle Physics, Santa Cruz, CA 95064 
(Dated: September 3, 2010) 

Abstract 

A generic prediction in the paradigm of weakly interacting dark matter is the production of rela- 
tivistic particles from dark matter pair-annihilation in regions of high dark matter density. Ultra- 
relativistic electrons and positrons produced in the center of the Galaxy by dark matter annihilation 
should produce a diffuse synchrotron emission. While the spectral shape of the synchrotron dark 
matter haze depends on the particle model (and secondarily on the galactic magnetic fields), the 
morphology of the haze depends primarily on (1) the dark matter density distribution, (2) the 
galactic magnetic field morphology, and (3) the diffusion model for high-energy cosmic-ray leptons. 
Interestingly, an unidentified excess of microwave radiation with characteristics similar to those 
predicted by dark matter models has been claimed to exist near the galactic center region in the 
data reported by the WMAP satellite, and dubbed the "WMAP haze" . In this study, we carry out 
a self-consistent treatment of the variables enumerated above, enforcing constraints from the avail- 
able data on cosmic rays, radio surveys and diffuse gamma rays. We outline and make predictions 
for the general morphology and spectral features of a "dark matter haze" and we compare them to 
the WMAP haze data. We also characterize and study the spectrum and spatial distribution of the 
inverse Compton emission resulting from the same population of energetic electrons and positrons. 
We point out that the spectrum and morphology of the radio emission at different frequencies is 
a powerful diagnostics to test whether a galactic synchrotron haze indeed originates from dark 
matter annihilation. 
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I. INTRODUCTION 



A compelling paradigm for the particle nature of the dark matter is that of Weakly 
Interacting Massive Particles, or WIMPs 

HQ 

. Although the Standard Model of particle 
physics does not encompass a viable particle dark matter candidate, WIMPs are predicted 
to exist in several well motivated extensions. These include weak-scale supersymnietry jsj. 



). WIMPs 



models with universal extra dimensions [4] , and many others (for reviews see 
are massive particles with masses near the electro-weak scale, and are typically charged 
under weak interactions. General arguments indicate that WIMPs in thermal equilibrium 
in the very early universe would freeze-out, decoupling from the thermal bath when the 
temperature dropped below a fraction (typically ~l/20) of their mass [ij]. The remaining 
dark matter particles would then populate the universe with a relic density which is of the 



same order as dark matter on cosmo logical scales [5|-|8|. This 'WIMP miracle" jo] warrants 
extensive investigation, due to the possibility of observing the particle debris stemming from 
the occasional dark-matter pair-annihilation event in today's cold universe. 

The probability of two dark matter particles (x) pair-annihilating into observable stan- 
dard model particles is proportional to the thermally averaged pair-annihilation cross section 
times the relative particle velocity, (av) , multiplied by the local particle dark matter number 
density squared, n^. The first quantity, (crv), can be inferred from the requirement of having 
a relic abundance ^ 1/ (av) on the same order as the universal dark matter density, i.e. 
^DM ~ 0.24 [10]. The latter quantity is given by = (pdm/''^x)^' "where indicates the 
mass of the dark matter particle x- Most dark matter annihilation processes are therefore 
predicted to occur in regions with a large dark matter density. Both intuition and detailed 



results from N-body simulations (see e.g. Ill-ll3l|) indicate that the center of the Milky Way 



galaxy is likely the brightest local dark matter annihilation site (barring the possibility of 
highly concentrated local dark matter clumps [3]). Among the possible signatures of dark 



matter annihilation, extensive studies have focused on gamma rays (see e.g. 15^, uM and 

species which have 

the benefit of carrying directional and spectral information. 

Other stable standard model particles produced in the pair annihilation of dark matter 
include charged cosmic rays such as electrons and positrons (e^) as well as (anti-)protons 
and (anti-)deuterons. These charged species scatter off of magnetic field irregularities in the 



Galaxy, losing energy and diffusing before reaching the Earth. The random-walk propagation 
of charged cosmic rays in the Galaxy is usually described with a diffusion-loss equation and 
solved numerically [18] or semi-analytically |l9l . |20|. A pedagogical review of cosmic-ray 



propagation and interactions in the Galaxy is given in Ref. 



2l|. 



The possibility of detecting an anomalous spectral feature in the flux of leptons, which 
could be traced back to the pair-annihilation of particle dark matter, has been long discussed 
(for early studies see e.g. Ref. 22|-|25|). This scenario has recently gained great momen- 
tum after results from the Pamela space-based antimatter detector reported an excess of 
high-energy (10-100 GeV) positrons over the assumed background of secondary positron 



production by inelastic cosmic-ray interactions |26|. Tantalizingly, an excess (namely an 
anomalous "bump" ) in the flux of electrons plus positrons was also recently reported by the 
balloon-borne experiments ATIC [27! and PPB-BETS [28] at energies in the range of several 
hundred GeV. This excess was subsequently not confirmed by data from the Fermi Large 
Area Telescope (LAT), which, however, did indicate a much harder spectrum for the flux 



at high energy than previously assumed 



29 



30| . This implies that the positron deficit re- 



ported by Pamela is at an even greater contrast with the standard expectation from cosmic 
ray models. The Fermi results agree with the low-energy range of other determinations of 
the flux [31]. 

Although astrophysical sources such as pulsars 
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3214361] and supernova remnants 



371, 



have been shown to provide a possible explanation for the positron fraction anomaly, 
dark matter models have been formulated which can fit both cosmic ray data and account 
for the positron excess (for a list of references see e.g. 3^)- In general, dark matter 
models that account for the Pamela anomaly need to have a large pair-annihilation cross 
section compared to the estimates from the standard thermal-relic calculation. Furthermore, 
a mechanism or conservation law must be invoked to ensure that excess antiprotons are 
not produced at a detectable level, as the Pamela data place stringent constraints on any 



antiproton excess [39[. Dark matter models which match both the Pamela and Fermi-LAT 
data commonly feature pair-annihilation modes dominated by positron production, and 
have rather large masses (~ 1 TeV). Recently, Ref. |40j] also pointed out that a wino-like 
candidate with a somewhat lighter mass can also account for the noted cosmic ray anomalies. 
We will entertain this possibility for one of the particle dark matter models considered below 
(although recent data on gamma-ray observations from local dwarf galaxies put stringent 



constraints on this category of models (4l|). Winos dominantly pair-annihilate into W^W~ 
pairs, with a significant subsequent production of energetic e^. 

While the source term for the injected by dark matter pair annihilation traces the 
dark matter density profile squared, the resulting non-thermal population of electrons and 
positrons settles to an equilibrium configuration only after losing energy and spatially dif- 
fusing through scattering off of both the regular and turbulent components of the galactic 
magnetic fields. At high energy, the dominant energy loss processes are synchrotron and 
inverse Compton scattering, which create secondary radiation at radio frequencies as well 
as X-ray to gamma-ray energies, respectively (for a recent review of multi-wavelength emis- 
sions from dark matter annihilation see Ref. j^]). This secondary emission is a common 
denominator to any WIMP model. In particular, dark matter models that also explain the 
Pamela positron excess generically give very large signals, which are typically barely com- 



patible with the extragalactic gamma-ray background, as shown e.g. in 43 



. Constraints 



on dark matter models from the secondary radio and X-ray emission from produced in 
dark matter annihilation in extragalactic objects such as dwarfs and galaxy clusters has 



been studied in a number of recent analyses, see e.g. 



46|-|52|. 



A. The role of synchrotron emission 



WIMP annihilation into charged leptons should produce synchrotron emission throughout 
a large spherical halo around the Milky Way. Placing the injection energy on the weak scale 
should create synchrotron emission at radio frequencies, since t^ galactic magnetic fields 
are of C ~ 1 — lOyuG. This fact was envisioned long ago, see e.g 



53 



56|, and, in conjunction 



with available radio data, was recently used to put severe constraints on the detectability of 
a gamma-ray signal from the galactic center jST] as well as on the viability of dark matter 
models that could explain the Pamela excess [58|. 

Interestingly, Finkbeiner pointed out in 2004 that a residual, unaccounted-for ra- 
dio "haze" could actually be present in the data reported by the Wilkinson Microwave 
Anisotropy Probe (WMAP) 59|], and that this radio emission could in principle be ex- 
plained precisely with the radio haze predicted in typical WIMP modelsO This point was 



^ Throughout this paper we refer to the observed WMAP signal as the "WMAP Haze" and our simulated 
results as the "dark matter haze" or "synchrotron haze" 
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7Q|, 



further elaborated upon in Ref. |60l-l66l|. which essentially confirmed that the morphology 
and intensity of the haze from WIMP annihilation might match the observed radio excess 
in the WMAP data. Other studies have shown that pulsars can create an additional lepton 
component which adequately traces the morphology of the WMAP haze [g?!. We note, 
however, that recent analyses by the WMAP team finds no evidence for a haze emission 
from the galactic center region in the polarization data [68j]. Furthermore, the extraction 
of a haze may depend sensitively on the assumption that the 408 Mhz synchrotron map is 

ous 



morphologically similar to skymaps in the WMAP band [69!]. Due to both the anoma 
positron results, and in view of the recent successful launch of the Planck satellite 
there is considerable interest in models with a large annihilation rate into leptons. Thus, we 
consider this a timely moment to critically re-assess the diffuse radio emission from WIMP 
annihilation in the central regions of the Galaxy and to outline possible diagnostics that 
will allow to disentangle it from other possible astrophysical sources. 

The secondary radio emission from dark matter annihilation depends on both the injec- 
tion spectrum and the pair-annihilation rate of dark matter - quantities which are set by 
specifying the particle dark matter model at hand. In addition, the intensity and spatial 
distribution of the WIMP synchrotron halo depend on three astrophysical quantities: (1) 
the dark matter density distribution, (2) the diffusion setup, and (3) the magnetic field 
setup. While information on (1) can be gained from N-body simulations, with educated 
estimates for the role of baryons and stars in shaping the dark matter density distribution, 
a self-consistent guess for (2) and (3) is a difficult task. Thus, the production of accurate 
predictions involves not only factoring in as many astronomical and cosmic-ray data as pos- 
sible, but also treating with accuracy the distribution of astrophysical cosmic-ray sources 
and the diffusion and energy losses of galactic cosmic rays. 

The importance of the diffusion setup, including in particular the geometry of the diffusion 
zone, the magnitude of the diffusion coefficient and its rigidity dependence, as well as the 
spatial distribution of Galactic magnetic fields, cannot be understated in the context of 
calculating the radio emission from dark matter annihilation. For instance, the height of the 
diffusion zone determines at which galactic latitude the population responsible for the 
radio emission leaks out and is artificially truncated by the boundary conditions imposed to 
the diffusion-loss equation; similarly, the spatial distribution of magnetic fields determines 
the intensity of the synchrotron emission, which scales quadratically with the magnitude of 



the magnetic fields. 

In the present work we study the morphology of the synchrotron halo produced by elec- 
trons and positrons created via dark matter annihilation in the context of self- consistent 
diffusion models. In particular, we interface the most complete publicly available cosmic-ray 
propagation code, Galprop 18|, with the spectra from dark matter annihilation cal- 
culated by the DarkSUSY package j71|. Testing various WIMP annihilation pathways, we 
compare our observations to constraints from cosmic-ray observations (e.g. Boron to Carbon 
and ^°Be to ^Be). Imposing this set of observational constraints on the diffusion setup forces 
the morphology of the dark matter "haze" to specific configurations, which we outline in 
detail and compare to the WMAP data residual as analyzed in Ref . [72| . 

Finally, we explore the diffuse inverse Compton emission from the up-scattering of the 
inter-stellar radiation field, which stems from the same population of dark-matter-induced 

that would produce the radio haze. We then compare the predicted emission with the 
Fermi diffuse gamma-ray data, and note the ability for observations of 7-ray emission to put 
strong constraints on models of the WMAP haze. 

The outline of this study is as follows: in the next section we describe the "baseline" 
setup adopted for the benchmark particle dark matter models (Sec. Ill A|) and for the set 
of galactic dark matter density profiles (Sec. IIIB|) . In addition, we provide details on the 
diffusion setup for cosmic ray propagation we employ in this analysis (Sec. Ill C|) . We study 
how the morphology of the synchrotron emission from dark matter annihilation depends 
on the diffusion setup in Sec. HVl on the dark matter density profile in Sec. El and on the 
spatial distribution of magnetic fields in Sec. I VI I We then study the impact of matching the 
adopted diffusion scenarios to cosmic ray data in Sec. I VI It and outline the expected spectrum 
of radio emission and inverse Compton emission in the context of the self-consistent scenarios 
in Sec. I Villi and IIXI respectively. Finally, we give an overview and a discussion of our results 
and predictions in Sec. |Xl 



II. BASELINE SETUP 

In this section we define a basehne setup for the present study. In particular, we discuss 
the benchmark particle dark matter models we employ (sec. Ill Al) . and the set of Galactic 
dark matter density distribution profiles (sec. IIIBI) . Both the particle setups and the density 
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distributions are chosen with the intention of illustrating extreme cases, and of capturing 
the full range of possible outcomes for our subsequent calculations. We then discuss our 
treatment of diffusion and energy losses, specifying the details of the setup we employ here 
as well as the assumptions on the magnetic field spatial distribution in sec. Ill CI 

A. Dark Matter Particle Models 

The injection spectrum of electrons and positrons from particle dark matter annihilation 
depends on (a) the dark matter mass, which sets the energy scale for the spectrum, and (b) 
the relative probability for the dark matter to pair-annihilate into given Standard Model 
final states. These branching ratios determine the spectral shape of the injected e^, which 
is calculated from the spectral distribution of electrons and positrons resulting from a given 
dark matter particle state. These spectral distributions are calculated by interpolating 
results from Monte ^arlo simulations, a task which is carried out automatically in the 
DarkSUSY package {zil. Finally, the overall normalization is directly proportional to the 
thermally averaged pair annihilation cross section times the relative velocity, (crt;), for v 0. 

Electrons and positrons are produced in dark matter annihilation events with spectral 
features which depend on the production mechanism. In particular, they can be produced 
(1) promptly, (and with an energy Ef.± ~ m^) in the process XX ~^ e~^e~, (2) with a hard 
spectrum, if they originate from muon or tau lepton decays, or (3) with a soft spectrum if 
they stem from the decay of charged pions resulting from hadronization chains of strongly 
interacting particles or from the hadronic decay modes of non-strongly interacting particles 
(tau leptons, gauge and Higgs bosons). Several theoretically motivated dark matter particle 
setups predict dominant pair-annihilation into a pair of charged SU(2) gauge bosons (this 
is the case, for instance, in anomaly mediated supersymmetry breaking (tSI]). In this case, 
the predicted injection spectrum is a combination of the above mentioned soft and hard 
spectral shapes. The soft channels correspond to the hadronic decay modes of the W, 
and the hard channels to the W Ivi leptonic modes, which produce energetic either 
promptly or via lepton decay. 

In the present study, we try to encompass the largest possible variety of WIMP dark 
matter injection spectra, with a preference toward models that are theoretically and/or 
phenomenologically motivated. For consistency, we set the global cross-section for all dark 
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matter particles to be (erf) = 3 x 10~^^ cm^/s and evaluate the necessary boosts from this 
value in order for our models to match the intensity of the WMAP haze. We note that 
this implies large boost factors for the case of very massive dark matter particles to get a 
significant signal. 

Specifically, we consider a "vanilla" dark matter model which corresponds to realizations 
of minimal supersymmetric extensions of the Standard Model as a benchmark for a soft 

injection spectrum. To populate the halo with relatively low-energy from WIMP 
annihilation, we set the dark matter mass to 40 GeV. This corresponds to the lightest 
bino-like neutralino in models with gaugino unification at the grand unification scale. It is 
customary to assume, for light enough dark matter particles, a bb final state, given that in the 
case of supersymmetry, if both the ti channel is kinematically closed and decays into gauge 
or Higgs bosons are suppressed by the neutralino composition, helicity suppression dictates 
that the dominant fermion-antifermion final state corresponds to the heaviest fermion. We 
dub this model as Soft, from its spectral shape and mass range. 

We choose our opposite extremum benchmark to have a hard spectrum and a large 
mass. Specifically, we select a model with exclusive pair-annihilation into muons, i.e. XX ~^ 
fi~^fj^~. These models are motivated by successful fits to the Pamela positron anomaly data as 
well as to the Fermi-LAT data, as shown for instance in Ref . 7^ . A further requirement 
for this type of phenomenologically motivated setup is a rather large mass scale, which we 
fix here at 1.5 TeV. Theoretical frameworks where particle models with similar injection 
spectra, and comparably large pair annihilation cross sections were proposed for instance in 



Ref. 



and 



76| . We indicate this model as Hard because of its spectral appearance. 

Finally, to interpolate between the mass ranges and the spectral shapes, and to pick 
an alternative particle dark matter setup that explains the anomalous Pamela positron 



Model ID 


Mass (GeV) 


(o"u)/(cm^/s) 


Ann. Final State 


Soft 


40 


3 X 10-26 


bb 


Hard 


1500 


3 X 10^26 




Wino 


200 


3 X 10-26 


w+w- 



TABLE I. Physical properties of the benchmark particle dark matter models we consider in this 
study. 
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FIG. 1. Left: The spectrum (multiplied by energy squared) of electrons (and positrons) produced 
by tfie pair annifiilation of dark matter for our tliree bencfimarli models: Soft (featuring a ligtitest 
neutralino mass ofm^ = 40 GeV and a dominant annifiilation final state into bb), Wino (m^ = 200 
GeV, W~^W~ ) and Hard (m^. = 1.5 TeV, fi^fj,^). Rigtit: tlie same, but for ttie spectral energy 
distribution of gamma rays promptly produced at dark matter anniliilation, for tfie tliree benchmark 
models. 



excess (see however the recent results of Ref. 4l|]), we consider a model along the lines 
of those discussed in Ref. ^|. Specifically, we consider a mass of 200 GeV and model a 
wino-like neutralino dark matter as a particle with a 100% pair-annihilation branching ratio 
into W~^W~ pairs. We dub this model as Wino. Notice that we do not assume a pair- 
annihilation cross section at the level expected for wino-like dark matter [4o|, in order to 
compare directly the boost factors among the three models under consideration. The "boost 
factor" corresponding to a 200 GeV wino would be on the order of 10^. 

In Tab. [I] we summarize the physical properties of the benchmark dark matter models 
we employ in the present study: namely the particle mass, its thermally averaged pair- 
annihilation cross section times velocity and its dominant annihilation final state. Fig. [H 
left, shows the differential injection spectrum (prior to propagation and energy losses) for 
the three models under consideration, where the flux is multiplied by energy squared and in 
units of GeV per annihilation. The spectral properties discussed above are clearly visible in 
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the plot, including: the soft e spectrum from pion decay in the case of the bb final state, the 
mixed soft-hard spectrum of winos, with a peak corresponding to the leptonic production 
modes, and the ~ spectrum from fi^fi~ decay. 

We also show in the right panel of Fig. [1] the prompt gamma-ray spectrum produced in 
the same annihilation event that produces the e^. While the bb and the W^W~ spectra are 
dominated by the two-photon decay of neutral pions, the spectrum of gamma rays from the 
muon final state entirely stems from final state radiation, and scales as (see e.g. \(\)- 
This "primary" emission will also play an important role when estimating the overall gamma- 
ray emission from dark matter annihilation in the central regions of the Galaxy; although 
the spectral details and the normalization of this component are not directly related to the 
secondary emission which is the core of the present study. 



B. Dark Matter Galactic Density Profiles 

The intensity of indirect dark matter signals from both prompt and secondary emission 
(inverse Compton and bremsstrahlung at high energy, synchrotron at low energy) are highly 
dependent on the distribution of dark matter density throughout the Galaxy. In the case of 
annihilating dark matter, the total interaction rate goes as the integral over the local dark 
matter density squared. We note that in the case of decaying dark matter, the annihilation 
rate goes instead as the integral over the local density. 

In this study we restrict our evaluation to several spherically symmetric dark matter 
profiles. We fix the local dark matter density at the solar position {Rq ~ 8.5 kpc) to be 
equal to p{Rq) = 0.389 GeV/cm^ [78:|. We express each of the dark matter density profiles 
as functions f{x) of the dimensionless variable x = R/a, where R is the distance from the 
center of the Galaxy and a is a scale radius for the density profile. 



Motivated by recent high resolution collisionless N-body simulations j79|, we adopt the 



best-fitting Einasto density profile obtained from the results of the Aquarius simulation 



2 

/Einasto(a;) = exp (x" - 1) , a ~ 0.17 (1) 



a 



and a profile extrapolated from the Via Lactea II simulation 12 



/viaLactea(a;) =X-^(x + l)-=^+^ 7 ~ 1.24. (2) 
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The scale radius for /Rinasto is assumed to be aAquarius — 22 kpc, a value which a' 



what found in Ref. 
Ref. 



80[, while that for /via Lactea is set to ayia Lactea ~ 28.1 kpc 



so matches 



points out that their results can be fitted with the prototypical Navarro- Frenk- 



White density profile 



8l| . for which: 

/NFw(a^) 



X 



-Hx + l) 



(3) 



and where we assume anew — 22 kpc. To illustrate the effect of the scale radius, we will 
also consider values for anew as large as 100 kpc and as small as 10 kpc. 

Physical processes related to the accretion of baryons clumped within the central regions 
of the Galaxy 
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can lead to a smoothing of the central cusps of the profile obtained 
in coUisionless dark matter simulations, and to the formation of large cores. The resulting 
density profiles can be modeled with profiles similar to those proposed in Ref. 8^, 85 1, which 
we indicate as a Burkert profile, with a functional form: 



/Burkcrt(x) = (1 + x)-\l + X^' (4) 

Galactic dynamics data point to OBurkert — 11-6 kpc. For illustrative purposes, we will 
consider a range of scale-sizes, with aeurkert — 5 kpc, OBurkert — 20 kpc, and OBurkert — 50 kpc. 

Several models of the galactic dark matter distribution predict cusps of high dark matter 
density near the galactic center. In this study, we do not include such cusps. In Section 
V, we shall point out that the impact of dark matter cusps on the morphology of the dark 
matter synchrotron haze can be understood and extrapolated off of our non-cusped models. 



C. Diffusion Models 



We utilize the state-of-the-art cosmic ray propagation code Galprop |18l. l86l|. which was 
created to solve cosmic-ray transport in a self-consistent manner. Galprop contains modules 
to compute secondary to primary ratios, synchrotron radiation, inverse Compton scattering, 
vr'^-decay, and bremsstrahlung radiation, from the same astrophysical inputs. These calcula- 
tions include effects such as nuclear spallation and energy loss due to interactions between 
cosmic rays and the Galaxy's gas, radiation, and magnetic fields. 

Utilizing a 2D simulation, we choose a standard diffusion model (GALDEF 50p_599278)cl, 



http://galprop.stanford.edu/GALDEF/galdef_50p_599278 
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based on the latest public version of the Galprop code. We note that the parameter space 
for propagation of cosmic rays is very large, and thus we take default assumptions consistent 
with the above model unless otherwise noted. In Table III CI we show the default values for 
the parameters most relevant to our study. We will evaluate first order changes in each 
of these parameters, while keeping the remaining terms constant to the values indicated 
in the Table. We place upper and lower bounds on each parameter, the values of which 
are not necessarily physically motivated. Instead, these parameters will serve to illustrate 
the effect of each individual parameter on the diffusion of high energy lepton sources. In 
Sec. IVIIt we will examine examples of the interdependence between diffusion parameters, as 
well as limiting constraints on possible model parameters. In order to model the galactic 
ICS emission from our lepton population (see Section IX), we employ the default Galprop 
model for the interstellar radiation field (ISRF) 87 1. 

We note that the assumed galactic diffusion setup may not be constant throughout the 
galaxy. Specifically, diffusion may be significantly different near the galactic center, due, 
for instance, to high gas and dust densities, which may affect the morphology of both 



synchrotron signals 



89|. In Section VII we 



as well as local cosmic ray abundances 
discuss the effect of these parameters on our simulated synchrotron morphologies. 

In order to determine the e^ fluxes from each of our dark matter annihilation pathways, 
we have modified the Galprop code to accept the input spectrum directly from the 
output of the DarkSUSY code [90|. The DarkSUSY code has been created to calculate the 
decay branching ratios and energy spectra for neutralino dark matter models as an input to 
the determination of dark matter direct and indirect detection rates. DarkSUSY is still the 
state-of-the-art tool in calculating the injection spectra of stable Standard Model particles 
which result from given particle dark matter setups. The Galprop code then interpolates 
between the input injection energy spectrum of e^, and calculates the approximate resulting 
flux for each necessary lepton energy. 



III. THE HAZE MORPHOLOGY: DEFAULT MODEL 

From our Galprop simulations, we calculate the spatial distribution of synchrotron ra- 
diation over the line of sight from the solar position. In order to determine the intensity 
of synchrotron radiation which stems from products of dark matter annihilation, we create 
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Parameter 


Default Value 


Lower Bound 


Upper Bound 


Alternative 


XXdlilUo 










Height 


4 kpc 


2 kpc 


16 kpc 




Dn 


5.0 X 10^^ cm^s"^ 


1.0 X 10^^ cm^s"^ 


1.0 X 10'^° cm^s~^ 




galaxy 


11.69e-^/i°'^P'=e-^/2kpc^G 


11.69e-''/i°'^P^e-^/i'^PVG 


10/x G 


11.69e-''/^°''P'^e-^/8''P^/iG 


Vq 


25 km s""*^ 


km 


100 km 




Convection 


Disabled 




100 km 


50 km 



TABLE IL List of important Galprop parameters used throughout this paper. Our default values 



are listed in the 2nd column, while important parameter space checks employed throughout this 
paper are listed in subsequent columns. 




10 20 

Galactic Latitude (Degrees off Center) 

FIG. 2. Simulated dark matter haze for three dark matter channels (4-0 GeV XX -> ('Soft, black 
solid), 200 GeV XX -> W+W' fWino, red dotted), 1500 GeV xX -> f^^fJ'^ (Hard, blue dashed)) 
plotted against the WMAP haze as determined by Dobler and Finkbeiner [72]. 

simulations both with and without dark matter for each set of diffusion parameters. Sub- 
tracting the dark matter signal, pixel by pixel, from the full simulations leaves a residual 
which stems purely from dark matter annihilations. 

For the observed WMAP haze, we use the fits derived from CMB5 in Dobler and 
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Finkbeiner 72|]. However, we fit the foregrounds over tlie entire sky, so that the analog 
of Eq. (18) in Dobler and Finkbeiner 72| becomes: 



rn = rpss + o-hH (5) 



The resulting map is then mean-subtracted to remove any unwanted residuals stem- 
ming from inadequate subtraction of the extragalactic isotropic background, and the data is 
displayed in unmasked regions (corresponding to the Dobler and Finkbeiner 72| mask) av- 
eraging the longitudinal values over the angular window -10° < 1 < 10°, in order to decrease 
any spurious noise effect. 

In Figure [21 we show the synchrotron overabundance at 23 Ghz attributable to dark 
matter in our default simulation, compared against observations of the WMAP haze. We 
multiply the output from our dark matter model by a fixed constant in order to normalize 
the diffuse radio emission ("haze") to the WMAP haze at a position of 10° latitude. This 
normalization accounts non-discriminantly for any change in the dark matter annihilation 
rate, but does not model the spatial characteristics which may be attributable to e.g. Som- 
merfeld enhancement or dark matter substructure. Since the boost factor is a well known 
uncertainty in calculating dark matter annihilation rates, we only concern ourselves here 
with matching the latitudinal distribution of dark matter, which depends primarily on the 
much better constrained global distribution of dark matter. 

In examining Figure [2l we note a significant contrast between the spatial distribution 
of the observed WMAP haze and our default simulation parameters. While the WMAP 
haze falls off by approximately a factor of 6.7 between 7° and 30° of latitude, our default 
parameters show a much steeper trend, falling off by factors of ~61 for our Soft model, 
~31 for our Wino model, and ~21 for our Hard model. Because this discrepancy does 
not depend on the overall annihilation rate, it is difficult to explain with variations in the 
dark matter annihilation setup. Instead, this discrepancy would most likely be explained 
by astrophysical parameters, and thus we first examine the role of our chosen diffusion 
parameters in determining the spatial distribution of the dark matter haze. 
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FIG. 3. Simulated dark matter haze for four dark matter density profiles: NFW (black line), Via 
Lactea (red dotted), Einasto (blue dashed), and Burkert (green long dashed) for models in which 
no diffusion occurs and the secondary electrons stay exactly at the point of annihilation. The 
parameters for each dark matter profile are given in Sec. 177-Bl 

IV. THE HAZE MORPHOLOGY: ROLE OF DIFFUSION PARAMETERS 

In order to illustrate the importance of our diffusion setup in determining the morphology 
of the WMAP haze, we show in Figure [3] the simulated WMAP haze for models in which 
we have entirely disabled diffusion (setting the diffusion constant, Alfven and convection 
velocity to 0). In this model, the dark matter synchrotron overabundance is a tracer of the 
position of dark matter annihilation and is identical to both the integral over the line of 
sight of the dark matter density squared, and the integral over the line of sight of prompt 
emission from dark matter annihilation. Interestingly, we note that in this case, the ratio 
between 7° and 30° is actually closer to unity than in our default simulation. This suggests 
that cosmic ray diffusion does not help equilibrate medium and high latitude emission, but 
instead greatly increases secondary emission near 10° by moving a large number of electrons 
from dark matter near the galactic center up to small galactic latitudes. 

Since changes in the diffusion setup can produce counter-intuitive changes like those 
shown above, it is important to create models of the dark matter haze over wide ranges of 
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FIG. 4. Simulated dark matter haze for three dark matter channels (4-0 GeV XX "> ('Soft, left), 
200 GeV XX -> W'^ W'^ (Wino, center), 1500 GeV xX -> A^^A*^ (Hard, right)) for illustrative 
choices of four important diffusion parameters (diffusion constant (top), diffusion region height 
(2nd row), Alfven velocity (3rd row), and convection velocity (bottom) plotted against the WMAP 
haze as determined by Dobler and Finkbeiner ]72l ]. In each plot, our default simulation is shown 
in solid black, with changes in each parameter shown in dotted red and dashed blue. 
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diffusion setups, in order to determine whether any of these provides a reasonable match 
to the WMAP haze. To determine the first order effects of changes in the Galprop diffu- 
sion setup on the spatial distribution of the dark matter haze, we modulate each diffusion 
parameter independently around its default value while keeping the remaining diffusion pa- 
rameters constant. In Figure IU we plot the synchrotron overabundance attributable to dark 
matter at 23 Ghz for changes in the diffusion constant (Dq), diffusion region height [Z^ax), 
Alfven velocity (vq), and convection velocity {vconv) for each of our dark matter annihilation 
pathways. We again normalize each model to match the WMAP haze at 10° latitude. In 
each plot, the black solid line shows our default propagation setup, while the various colored 
lines show illustrative changes in a single propagation parameter. 

As shown in Figure [3l cosmic ray diffusion actually increases the falloff in emission be- 
tween 7° and 30°. However, on the opposite end of the spectrum, we would expect completely 
isotropic emission in the limit of infinitely large diffusion. Interestingly, we find that the 
standard diffusion constant of Dq = 5.0 x 10^^ cm^s~^ lies near a maximum in the ratio of 
interior (i.e. low-latitude) synchrotron production to high latitude emissivity. Thus, moving 
away from this value of the diffusion constant in either direction helps bring around a better 
match to the WMAP haze morphology. We find the effect of limiting diffusion from the 
galactic center to low latitude regions to be dominant, as low values of the diffusion con- 
stant tend to greatly increase the ratio of high latitude to low latitude synchrotron radiation, 
providing a best match to the WMAP haze for our Soft model at Dq ~ 1.0 x 10^^ cm^s~^ 



consistent with the results of Hooper et al. 6l|), while our Wino and Hard models require 



even lower diffusion constants. High values of the diffusion constant provide a slightly better 
match to the WMAP haze, but the effect is never strong enough to provide a reasonable 
match, even for extremely large diffusion constants (Dq = 1.0 x 10^^ cm^s~^ or greater). 

Changes in the height of the diffusion zone (2nd row) have an interesting effect on dark 
matter models of the haze. For smaller diffusion zones {Zmax = 2 kpc), we see a rapid, 
and functionally different, decrease in the dark matter haze starting at approximately 10° 
latitude. This is due to the fact that from the solar position, a latitude of 10° implies a 
height of 1.47 kpc above the galactic center. Any observation at higher angles will quickly 
lose large quantities of emitting synchrotron radiation to the confines of the simulation 
zone. This restricts the normal method of making corresponding changes in the diffusion 
constant and simulation height (as is commonly done for studies of the local primary to 
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secondary ratios). However, from the simulation for "Lmax = 16 kpc, we find that even 
large increases in the diffusion zone beyond our default value create only small changes in 
the output haze. Furthermore, increases in the height of the simulation region will require 
corresponding increases in the diffusion constant. Thus, we find that a diffusion height of 
4 kpc is an optimal value for simulations of the WMAP haze up to 30° latitude. 

The Alfven velocity can affect the angular distribution of the synchrotron haze by altering 



the diffusion of dark matter secondaries in momentum space 9l|, l92| . The diffusion through 
momentum space varies as the square of the Alfven velocity: 

Dp, oc (6) 

and thus changes which decrease the Alfven velocity have little effect on the angular 
distribution of the haze. However, large Alfven velocities can create very effective diffusion 
through momentum space, which either propel particles up to high latitudes, or restrict 
particles from entering even the low latitude simulations regions. This mimics the improve- 
ments seen in diffusion constants that are both higher and lower than our standard setup. 
We find that an Alfven velocity of approximately 100 km s^^ can successfully match the 
morphology of the WMAP haze for each of our dark matter particle physics models. 

While diffusion is the dominant driving force for spatial diffusion of cosmic rays, galactic 
convection can also play a significant role. As we have convection disabled in our base model, 
we test the parameter space for models with stronger convection. In the bottom row, we test 
several convection velocities which linearly vary with the height above the galactic plane. We 
note that changes in the convection velocity do not greatly change the spatial distribution of 
the dark matter haze, and the variation they do create goes in the wrong direction to match 
the observed haze morphology. This is due to the fact that convection is normally strongest 
at higher latitudes, and thus efficiently moves high energy leptons out of our simulation. We 
also note that the effect of convection is strongest for soft dark matter models, as the effect 
of convection is most important between 1-10 GeV, which is a sweet spot for e^*" production 
from our soft model. 

Lastly, we note that a well known degeneracy exists between changes in the diffusion 
constant and the scale height of the diffusion region. While a change in the direction towards 
low diffusion constants and low scale heights are unlikely to create a match to the WMAP 
haze, as scale heights less than 4 kpc cut off emission above the galactic center at latitudes 
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FIG. 5. Simulated dark matter haze for three dark matter channels (40 GeV XX "> (Soft, left), 
200 GeV XX -> W'^ W'^ (Wino, center), 1500 GeV xX -> /^^/^^ (Hard, right)) for illustrative 
choices of the diffusion constant and simulation height which keep their ratio constant. In each 
plot, the black solid line shows our default propagation setup, while the various colored lines show 
our parameter space models. 




near 30°. However, in Figure HI we see that increases in both the diffusion constant and 
the height of the diffusion region increase the intensity of high latitude synchrotron relative 
to emission at 10° below the galactic center. This raises the possibility that corresponding 
changes in these parameters may create better dark matter fits to the WMAP haze, while 
preserving the local cosmic ray density. 

In Figure [5] we show the resulting dark matter haze for several parameter choices which 
leave the ratio of Dg/height constant. We find that larger diffusion constants and heights 
do increase the relative intensity of high latitude synchrotron radiation. However, we do 
not begin to see a good match until we employ extreme cases such as a diffusion height 
of (z=32 kpc), which itself still falls short of reproducing a match to the WMAP haze 
morphology, with reasonable matches occuring only in the cases of our harder spectra. We 
strongly caution the reader, however, about the accuracy of the synchrotron morphology 
for this extreme height, since the galactic magnetic field is almost completely unknown in 
this region. Nevertheless, as this is an important degeneracy condition in our work, we will 
examine it further in Section IVIII 
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V. THE HAZE MORPHOLOGY: ROLE OF DARK MATTER DENSITY PRO- 
FILES 



The distribution of dark matter secondaries should depend strongly on the galactic dis- 
tribution of dark matter annihilations. Throughout this paper, we model only the smooth 
component of dark matter, normalizing the result to the synchrotron haze in order to ac- 
count for uncertainties in the dark matter annihilation rate. Since we are concerned with 
annihilating dark matter, we note that the number of annihilations at any point in space 
varies with the square of the local dark matter density. To test the effect of various dark 
matter density profiles on our simulated haze, we examine several common dark matter dis- 
tributions in addition to our default NFW profile. We test the functional forms given by the 



Via Lactea simulation (R^c = 28.1 kpc) |12i]. an Einasto profile (R^c = 20 kpc, a = 0.17) 80|. 
and a Burkert profile (R^c = 11-6 kpc) SJ, l85|] . For all models we adopt a solar distance 



R© = 8.5 kpc and a local dark matter density of 0.389 GeV cm~^ 78|. In Figure El we have 
showed the normalized synchrotron intensity when the motion of dark matter secondaries is 
neglected. We note that this is identical (within a normalization constant) to the square of 
the dark matter density integrated over the line of sight for each of our assumed profiles. 

In Figure [6l we show the resulting dark matter haze for several choices of the NFW scale 
radius (top), various dark matter density profiles (middle), and several choices of the scale 
radius for the Burkert profile (bottom). As in previous plots, the results are normalized to 
the synchrotron haze at 10° latitude and at a frequency of 23 Ghz. Our default simulation 
is shown in black in the top two rows, but is not shown in the bottom row. 

We find that even large changes in the NFW profile scaling radius (top) have only a 
small effect on the morphology of the dark matter haze. For a scale radius of 100 kpc the 
ratio in synchrotron emission between 7° and 30° is still 37 for our Soft spectrum, 21 for 
the Wino spectrum, and 15 for the Hard spectrum (as opposed to 6.7 for the WMAP 
observation). This again illustrates the point made in Figures [3] and H] (top) - that the 
majority of synchrotron radiation seen at 30° is derived from leptons produced near the 
galactic center, which then diffuse to high latitudes. 

This analysis is again confirmed upon examination of several dark matter density profiles 
(middle). The NFW, Via Lactea, and Einasto profiles all predict central regions with high 
dark matter density, and all three profiles result in very similar morphologies for the dark 
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matter haze. This similarity comes despite the very differing DM densities near the galactic 
center in these profiles - with the average density throughout the inner 0.5 kpc of the Einasto 
profile being approximately twice as large as the Via Lactea profile and four times as large 
as the NFW profile. Since leptons produced in this region diffuse to high latitudes and 
dominate the dark matter synchrotron signal, these density profiles greatly affect the cross- 
section normalization in each model, but do not control it's morphology, which is governed 
instead by the diffusion setup. This allows us to infer that the addition of an additional even 
cuspier dark matter component into a profile with high central density will only significantly 
affect the assumed dark matter cross-section, with little effect, however, on the resultant 
synchrotron morphology. 

However, the Burkert profile, which does not include a dense galactic center, gives a 
spatial variation which is much more consistent with the WMAP haze. Because the Burkert 
profile makes dark matter annihilation at high latitudes relatively more important compared 
to annihilations at the galactic center, it creates a much flatter distribution similar to the 
results shown in Figure El However, we note that removing the dense galactic center region 
eliminates the vast majority of our dark matter induced synchrotron radiation. Thus, to 
match the intensity of the WMAP haze with the Burkert profile, we require much higher 
normalization constants (~3010 for the Soft dark matter model, ~2350 for the Wino model, 
and ~37700 for the Hard model) 

As a result, to create a match to the WMAP haze through changes in the dark matter 
density profile, one might invoke a distribution like the Burkert profile which increases the 
importance of annihilations occurring at high galactic latitudqj. In Row 3, we show the 
simulated dark matter haze for several values of the scale radius in the Burkert profile. 
Unlike the case of the NFW profile (top), the increased dependence on local dark matter 
annihilations in the Burkert profile creates a morphological dependence on the scale radius. 
For low scale radii like Ksc = 5 kpc, we see a sharper cutoff that does not match the 
morphology of the WMAP haze. In contrast, we see very good matches for scaling radii 
between 11.6-50 kpc, depending on the employed dark matter particle setup. 



^ Alternatively, we note that a galactic substructure distribution with a broad radial distribution of subhalos 
might also effectively provide a boost at large radii which can mimic a similar DM density distribution to 
the cored Burkert profile 
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FIG. 6. Simulated dark matter haze for three dark matter channels (40 GeV XX "> (Soft, left), 
200 GeV XX -> W'^ W'^ (Wino, center), 1500 GeV xX -> A*^/^^ (Hard, right)) for illustrative 
choices of the dark matter density profile (top) and scale radius in the NFW profile (bottom). In 
each plot, the default simulation is shown in solid black. 



VI. THE HAZE MORPHOLOGY: ROLE OF GALACTIC MAGNETIC FIELDS 



A final important parameter in determining the spatial distribution of synchrotron radi- 
ation is the morphology of galactic magnetic fields. Changes in the galactic magnetic field 
morphology may produce reasonable matches for the WMAP haze because magnetic fields 
(1) allow for variations in the spatial morphology of synchrotron radiation without greatly 
affecting cosmic ray or 7-ray observations, and (2) are rather poorly constrained away from 
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the solar position. However, we note that changes in the galactic magnetic field strength will 
not only affect the dark matter Haze, but will alter the intensity and morphology of all galac- 
tic synchrotron foregrounds, possibly eliminating or increasing the need for a WMAP haze. 
In this analysis, we artificially hold the observed WMAP haze constant through changes in 
the magnetic field intensity, as we only aim to show that known uncertainties in the galactic 
magnetic field could greatly alter the morphology of dark matter matter driven synchrotron 
emission. 

In Figure [71 we evaluate the dark matter haze for several models of the magnetic field 
which follow the functional form of 

5 = 5oexp(-r/ro -zjzQ), (7) 

for various values for zq. In addition, we model a fiat (technically: scale radius of 99.9 kpc in 
)oth z and r) component equivalent to that employed in previous studies, e.g. Hooper et al. 
6ol |. Again, the black line shows our base profile and we normalize all results to match the 
WMAP haze at 10° latitude. We find that employing a fiat magnetic field profile greatly 
increases the synchrotron radiation at high latitudes above the galactic center producing 
intensity ratios of 10.5 (Soft), 10.5 (Wino), 9.35 (Hard) between 7° and 30° in latitude. 
These values are more comparable to the ratio ~ 6.7 in the WMAP haze. However, we see 
that the fiat magnetic field profile does not match the sharp rise in the WMAP haze between 
7-10° in latitude and thus underproduces the haze in that region. The sharp magnetic field 
cutoff falls off very sharply at higher latitudes, while the smooth magnetic field model 
produces an intermediate case between the base and fiat profiles. We note that several 



previous attempts at dark matter models of the WMAP haze (e.g. 6l|) use a magnetic field 
which is spatially uniform, and show results that are consistent with those shown here. 

Since the spatial variation of the magnetic field far from the solar neighborhood is highly 
uncertain, all magnetic fields presented here are fairly reasonable choices. This creates a 
large uncertainty in putative dark matter models fitting the WMAP haze, as the morphology 
of an additional synchrotron component can be adjusted by at least a factor of two from low 
to high latitudes. In addition to the freedom of choosing a normalization constant to match 
the intensity of the WMAP haze, this allows a secondary normalization constant which 
effectively controls the morphology of the haze as a function of latitude. This fact can be 
interpreted in two ways: either as evidence that dark matter models can successfully describe 
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FIG. 7. Simulated dark matter haze for three dark matter channels (40 GeV XX "> (Soft, left), 
200 GeV XX -> W'^ W'^ (Wino, center), 1500 GeV xX -> /^^/^^ (Hard, right)) for illustrative 
choices of the galactic magnetic field. In this plot, the default simulation parameters are shown in 
black with steeper and smoother magnetic fields shown in the various colors and line types. 

an excess synchrotron haze, or that the existence of a dark matter model which matches the 
morphology of the WMAP haze is not convincing proof of the models correctness, because 
uncertainties in the magnetic field make the underlying diffuse morphology and spectrum 
difficult to pinpoint. 

VII. TOWARDS A SELF-CONSISTENT PICTURE: CONSTRAINTS FROM COS- 
MIC RAY DATA 

Until this point, we have simply made ad hoc changes to the diffusion, magnetic field, 
and dark matter density profiles without regard for constraints which may exist from other 
observations. In the remaining sections of this paper, we will much more closely examine the 
reasonability of our parameter space choices based on various other observable quantities. 
In particular, we outline the correlated ranges for the parameters influencing cosmic ray 
transport according to local constraints from cosmic ray data. 

As seen in Figure HI changes in both the diffusion constant and Alfven velocity can 
greatly affect the morphology of the dark matter haze, potentially allowing well-motivated 
dark matter particle physics setups to match the observed WMAP haze. However, changes 
in the diffusion setup do not affect only the leptons produced by dark matter annihilation, 
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FIG. 8. The Boron to carbon ratio (black solid / black triangles) and the ^^Be to ^Be ratio (blue 
dashed / blue squares) modeled at the solar position for three important cosmic ray propagation 
setups: Dq = 1 X cw? s"-*^ (l^ft), f^a = 100 km s'^ (middle), and our base model (right). 

but affect the morphology of all cosmic ray species as well. Since cosmic rays interact as 
they propagate in the Galaxy, ratios of secondary to primary cosmic ray species provide a 
valuable consistency check for the assumed cosmic ray diffusion models. We note that these 
same constraints do not exist for most changes in the dark matter profile or the magnetic 
field models, as neither have a large effect on the galactic primary to secondary ratios. 

In Figure [HI we show some of the simulated cosmic ray data (the Boron to Carbon Ratio 
and the ^°Be/^Be ratio) for the two diffusion models in Section IIV| which matched the 
morphology of the WMAP haze: (1) a low diffusion constant of Dq = 1 x 10^^ cm^ s~^ (left) 
and (2) a high Alfven velocity of Vq, = 100 km s~^ (middle)). We also show the nuclei fluxes 
for our default Galprop model (right). We note that the large shifts in diffusion parameters 
necessary to match the WMAP haze also have large effects on the simulated primary to 
secondary fluxes and ratios. Both models are definitively ruled out by current observations 
at levels much greater than lOcr. In our default Galprop model (which of course is refined 
in order to match observed cosmic ray data) we indeed find very close matches to current 
observations. 

However, in our analysis thus far, we have evaluated only first order changes in each 
diffusion parameter, while the rest of our diffusion setup is kept constant. It is entirely 
possible that some combination of changes will create a match between the dark matter Haze 
and WMAP haze, while maintaining the observed primary/secondary ratios. In Figure IH we 
show the two dimensional parameter space spanned by a reasonable range for the diffusion 
constant and for the Alfven velocity, which are two of the most significant parameters 
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affecting tfie dark matter liaze morpliology. For eacfi point in our parameter space, we 
calcufate a total from available cosmic ray ratio data 93l-l96| . The value of is calculated 
similarly to Cumberbatch et al. [97|, eq. 6], but without reducing the total variance: 



= S,Sf ^^^^^-^ (8) 

where D^, dij, and are the central values, la errors, and predicted abundances re- 
spectively of the ith datapoint with the jth experimental dataset. This value is represented 
in the z-axis of Figure M Contour lines trace out the regions where each ratio is consistent 
with the data to within three and five sigma. The tightest constraint comes from the Boron 
to Carbon ratio, which gives 3 a constraints limited to only a signal point in our gridding 
simulation, and 5 cr constraints on the diffusion constant above 4.0 x 10^® cm^s~^ and be- 
low 6.0 X 10^® cm^s~^. The other ratios also appear to have reasonable values around that 
same point, with a diffusion constant of Dq ~5.0 x 10^® cm^s~^ and an Alfven velocity of 
25 km s"^. 

In Figure ??, we show a similar two-dimensional plot of the diffusion constant compared 
the height of the diffusion region. We indeed note a direct correlation between higher dif- 
fusion constants and diffusion heights, which creates a best fit region which is significantly 
more extended than the two dimensional fit of the diffusion constant with Alfven velocity. 
We find no constraints from either the Ne or Mg ratios, and do not include these in the 
fit. However, reasonable simulation heights rule out diffusion constants larger than approxi- 
mately 1.5 X 10^^ cm^s~^. In Figure [5l we see that usage of a higher diffusion constant along 
with a simulation height of 8 kpc does lead to a reasonable improvement in the morphol- 
ogy of the dark matter haze, compared to WMAP observations, especially for our Hard 
annihilation spectra. However, even much larger simulation heights are unable to create a 
complete match to the haze, ruling out the possibility that this degeneracy condition alone 
is capable of explaining the morphological differences between the observed and simulated 
haze. 

We note that in this study we have only considered models employing diffusion parameters 
which are homogeneous throughout the diffusion zone. Physical intuition suggests that the 
high matter and radiation densities near the galactic center may locally decrease the diffusion 
constant and affect, for instance the re-acceleration of cosmic rays. While this may alter the 
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cosmic ray primary-to-secondary ratios, our models indicate it should have little effect on 
the assumed dark matter synchrotron morphology - as in the cases of the NFW, Einasto, 
and Via Lactea profiles, the dark matter induced lepton population is entirely dominated by 
annihilations occuring near the galactic center. As long as changes in the diffusion constant 
are not strong enough to change this assumption, then the specific details of diffusion near 
the galactic center are degenerate with assumptions on the dark matter cross section (for 
changes in the diffusion constant), the input spectra of dark matter produced leptons (for 
changes in the reacceleration parameters) and the dark matter density distribution in the 
galactic center region. In our models, the NFW, Einasto, and Via Lactea profiles all have 
very different central annihilation rates, and thus the morphological similarity between these 
models indicates that the synchrotron morphology is determined primarily by the diffusion 
setup outside of the galactic center, rather than by the details of the diffusion setup very 
near the galactic center. 

In summary, our findings do not rule out the possibility that changes to the cosmic-ray 
propagation parameters compatible with local cosmic-ray data allow for significant effects on 
the high-energy lepton equilibrium distribution from a baseline setup. However, we showed 
that drastic changes to the morphology of the galactic equilibrium number density are 
produced by parameter choices that tend to create substantial issues in the resulting ratios 
of primary-to-secondary cosmic-ray fluxes. 



VIII. FREQUENCY DEPENDENCE OF DARK MATTER HAZE 

One further test of the dark matter interpretation of the WMAP haze, originally proposed 
by Ando and Komatsu [oS^ , concerns examining the frequency dependence of simulated dark 
matter haze models. The main parameter space uncertainties for these models result from 
the dark matter annihilation pathway and the magnetic field model, which can produce either 
harder or softer spectra. In Figure fTTl we show the spectra of our default model, as well as a 
model employing a Burkert profile and the flat magnetic field model from Section IVIII The 
simulations are shown at 10° below the galactic center and normalized to match the observed 
data at 23 Ghz. The data shown here originates from all latitudes between 7-15 degrees, 
and averaged over the longitudinal strip -10° < / < 10° and are normalized at 23 Ghz to 
match the value at 10° below the galactic center. We show WMAP data at only 23 Ghz, 33 
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FIG. 9. Constraints from Nuclei fluxes and primary to secondary ratios as a function of Alfven 
velocity and diffusion constant. Darker areas indicate better fits to simulation data (lower x^)- 
The minimum falls around a diffusion constant of 5.0 x cm^s~^ and an Alfven velocity of 

25 km s^^ . Colored lines for each ratio show 3 and 5 sigma constraints on the dataset 



9S 



Ghz, and 41 Ghz, where the observed error bars do not dominate the residual signal. 
We point out two important results: 



1. the dominant contribution to the frequency dependence of the dark matter haze re- 
lates to the mass and annihilation pathway of the dark matter particle; Our three 
illustrative dark matter models show very different frequency dependences, which can 
help differentiate between each model, and 

2. the Wino model provides the best match for the spectral characteristics of the WMAP 
haze. 



We note, however, that we have not rigorously examined either the systematic or observa- 
tional errors in this measurement, and thus further study is needed to determine whether 
these frequency tests rigorously rule out any of the models presented here. Future measure- 
ments from the Planck satellite can be expected to provide an improved spectral analysis, 
making this technique and important multi-wavelength test of any future residual signal. 
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FIG. 10. Constraints from Nuclei fluxes and primary to secondary ratios as a function of the 
simulation height and diffusion constant. Darker areas indicate better fits to simulation data (lower 
X^)- A degeneracy exists along a line of constant diffusion constant divided by height. Colored lines 
for each ratio show 3 and 5 sigma constraints on the dataset. Some local maxima which are non- 
physical have been smoothed over. 

IX. COMPARISON TO INDIRECT DETECTION ON OTHER ENERGY SCALES 

While we have focused exclusively on the detectability of the synchrotron contribution 
of dark matter, we note that any additional primary lepton flux (be it from an annihilating 
dark matter candidate, pulsars, SN remnants etc.), is guaranteed to produce not only some 
synchrotron signal, but also signals at GeV energies due to inverse Compton scattering off of 
the interstellar radiation field of the Galaxy. We can thus put further constraints on our dark 
matter models of the WMAP haze by evaluating the expected 7-ray signature of each setup 
and enforcing constrains from recent Fermi data. We focus here on three morphologically 
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Flat D Field (10 /xG) 

Burkert Profile (R^ - 11 kpc) 
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a 

Synchrotron Frequency (Ghz) 

FIG. 11. Frequency dependence of the simulated dark matter haze for three dark matter decay 
channels (40 GeV XX -> bb (Soft, left), 200 GeV XX -> W+W' fWino, center), 1500 GeV XX 
-> fi^^L^ (Hard, right)) at 10 degrees below the galactic plane for our default and two models 
which provide reasonable WMAP haze matches: (1) base model (solid black), (2) Flat Magnetic 
Field (red dashed), (3) Burkert profile (blue long-dashed). The dotted region shows a range of 
WMAP frequency dependences for each degree in between 7-15 degrees below the galactic plane. All 
values have been normalized to match at 23 Ghz. 




interesting models: (1) our base setup (top), (2) a model with a flat magnetic fleld (mid- 
dle), and (3) a model employing a Burkert dark matter density proflle with Rc =11.6 kpc 
(bottom). In Figure [T2| we plot the expected Fermi signal due in each decay pathway due 
to both direct 7-ray production (green dot-dash), and inverse-Compton emission from scat- 
tering off of the galactic ISRF (red dash). For the total signal (black solid), we add emission 
from 7r° decay and bremsstrahlung emission from our default Galprop models, as well as 



actic isotropic background (Energy ) and point sources 



power law spectra for the extraga^ 
(Energy"^'^) (from inspection of [99j). 

We flnd that several of these models signiflcantly over-predict the observed Fermi 7-ray 
signal at intermediate latitudes. For our Soft and Wino models, we flnd that direct 7-ray 
production is the dominant DM signal at intermediate latitudes. The large boost factors 
necessary to recreate the intensity of the WMAP haze overproduce the observed Fermi 
signal in both our base and Burkert models, while the stronger magnetic fleld in our Flat 
B model allows the 7-ray emission to match observations. For our Wino models, we note 
an overproduction of direct 7-ray signals from our Burkert model, while the base and flat 
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magnetic field model both yield reasonable results. 

In the case of our Hard models, we find the dominant 7-ray production mechanism 
to stem from ICS of the ISRF. This signal overproduces the observed Fermi signal in our 
Burkert model, while it provides reasonable fits in our base and fiat magnetic field models. 

"or d ark matter annihilation in dwarf galax- 



100| (Figure 3, top) allows us to place a 



Stronger constraints also come from searches 
ies and clusters. The results of Abdo et al. 
cross-section limit of approximately 8 x 10^^^ cm'^ s^^ on bb annihilation spectra. This 
definitively rules out our Soft default model as it requires an annihilation cross-section of 
approximately 1.1 x 10~^^ cm'^ s~^. However, in the wino case (Figure 3, bottom), the best 
Fermi limits from dwarfs give a cross-section of approximately 1 x 10~^^ cm'^ s~^, which is 
well above our model Wino cross section of 8.3 x 10"^^ cm'^ s~^. In the case of our Hard 
spectrum, the best limits lie at approximately 5 x 10^^^ cm'^ s~^, which is well above our 
necessary cross-section of 1.2 x 10~^^ cm^ s~^. 

We also n ote s everal other relevant constraints on the ICS emission. Work by Papucci 



and Strumia 101| (Figure 3, NFW Profile) places limits of 4 x 10~^^ cm^ s~^ on our Soft 
model, 1.3 x 10^^^ cm^ s~^ on our Wino model, and 5 x 10^^^ cm^ s^^ on our Hard model - 
which gives better constraints on our wino model. Fermi models of the isotropic background 



102| (Figure 5, MSII-Subl limit) give 95% limits of 5 x 10 cm^ s ^ on our Soft model. 

As intuition suggests, the necessary cross-sections always decrease when we employ a 
stronger magnetic field model. In the case of our Soft annihilation spectrum, the flat 
magnetic field model lowers the necessary cross-section to just 1.83 x 10~^^ cm'^ s~^, which lies 
closer to the constraints imposed by Fermi observations. This yields the enticing possibility 
that low mass dark matter annihilation spectra with intensities able to match the haze will 
either be confirmed or ruled out by Fermi in the very near future. 

Notice that further constraints on dark matter annihilation and its emission from the 
galactic center region can be obtained by studying regions very close (on the order of a 
few arcsec) to Sag A*, under educated assumptions on the magnetic field intensity in that 



region 



57l |: in particular, synchrotron emission is severely constrained by both radio and 



57|] . Given the coarse spatial grid 



X-ray observations at small angular scales around Sag A* 
we employ in this study, and the fact that we are concerned with much larger angular scales, 
we do not discuss this type of constraints here. 
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FIG. 12. Expected Fermi Signal from each decay pathway for the base model (top), a flat B-field of 
10 nG (middle), and a model employing a Burkert profile with Rgc = 11.6 kpc bottom, and shown for 
three dark matter decay channels (4-0 GeV XX -> ('Soft, left), 200 GeV XX -> W^W" (Wino, 
center), 1500 GeV xx "> /^^^^ (Hard, right)). The red dash shows the expected Inverse Compton 
(ICS) emission based on the Galprop model for the galactic ISRF, while the green dot-dashed line 
is based on the direct j-ray production from a DarkSUSY model, and normalized based on the input 
dark matter density profile. The total emission (black solid) includes Galprop based emission for 
-decay and bremsstrahlung emission, as well as power law models for the extragalactic background 
(Energy'"^'^^ normalized to 1.1 x 10"^ MeV cm^^ s^^ sr^^ at 1 GeV) and point sources (Energy^'^'^ 
normalized to 3.1 x 10~^ MeV cm~'^ s^^ sr^^ at 1 GeV). The observed Fermi-LAT ^-ray flux is 
obtained from j^/. 
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X. DISCUSSION AND CONCLUSIONS 



WIMP annihilation in the Galaxy unavoidably yields high-energy electrons and positrons 
(barring the possibility of pair annihilations in neutrino pairs or particles outside the Stan- 
dard Model). In turn, this implies a diffuse galactic emission, whose morphology is affected 
by both the distribution of dark matter and charged cosmic-ray transport in the Galaxy. 
In this study, we outlined in detail how the form and nature of the WIMP radio halo de- 
pends upon the particle model for dark matter as well upon astrophysical parameters such 
as the galactic magnetic fields and the parameters entering the propagation of cosmic rays. 
We pointed out how cosmic-ray data very significantly constrain the morphology of the 
WIMP radio haze, and pointed out the importance of other secondary emission channels, in 
particular Inverse Compton. 

Given the interest in a possible excess emission at radio frequencies - the WMAP haze 
- we have shown that standard dark matter annihilation pathways and cosmic ray diffusion 
models do not provide a reasonable match for the morphology of the excess. Dark matter 
density profiles with high densities near the galactic center produce synchrotron morpholo- 
gies which are dominated by the diffusion of charged leptons away from the galactic center. 
However, changes in the diffusion setup do not allow one to reproduce the morphology of the 
WMAP haze while remaining consistent with local observations of cosmic ray primary to 
secondary ratios. Models which lack a dense dark matter distribution in the galactic center 
(e.g. the Burkert profile) require large cross-section enhancements which are ruled out by 
Fermi observations. 

We also found that the particle physics properties of dark matter models which seek to 
explain both the WMAP haze as well as higher energy observations such as those from 
the Fermi-LAT are beginning to be constrained. Low energy dark matter decay spectra 
such as our Soft model are unable to reproduce the intensity of the WMAP haze without 
overproducing the Fermi signal with 7-rays produced via direct annihilation. 

However many uncertainties still exist in the magnetic field morphology to allow viable 
matches for the observed excess in synchrotron emission. The magnetic field morphology 
primarily influences the synchrotron intensity and morphology, with very little effect on 7- 
ray or cosmic ray constraints. The best constraints on galactic magnetic fields would come 
from careful examination of the WMAP dataset which would determine the morphology of 
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any WMAP haze resulting from various magnetic field models. Such an analysis is beyond 



the scope o 
Finkbeiner 



this study, but has been examined in both Gold et al. [68| and Dobler and 
72| . with contrasting results. Furthermore, many parameters not mentioned 
in this work may be fine tuned to produce a reasonable match to all known constraints. 
However, the formidable size of this parameter space makes the existence of such matches 
unconvincing evidence that the WMAP haze arises from a dark matter source. 

Further inquiry in the nature of an excess diffuse radio emission i s cer t ainly warranted, 
especially since e.g. any additional IC signal in the high energy sky 1031 . Il04| implies the 
production of some microwave residual detectable in the microwave sky. Diffuse radio emis- 
sion, and constraints on it, can thus be regarded as a powerful multi- wavelength diagnostic 
of any claims of new anomalous signals. Lastly, future data collected by the Planck telescope 



70| will be imperative in determining the intensity of the WMAP haze, and in extending 



its spectrum to the higher energy regime. 
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